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Abstract 

Using event-by-event viscous fluid dynamics to evolve fluctuating initial density profiles from the 
Monte-Carlo Glauber model for U-fU collisions, we report a "knee'’-like structure in the elliptic flow as a 
function of collision centrality, located near 0.5% centrality as measured by the final charged multiplicity. 
This knee is due to the preferential selection of tip-on-tip collision geometries by a high-multiplicity 
trigger. Such a knee structure is not seen in the STAR data. This rules out the two-component MC- 
Glauber model for initial energy and entropy production. An enrichment of tip-tip configurations by 
triggering solely on high-multiplicity in the U-fU collisions thus does not work. On the other hand, using 
the Zero Degree Galorimeters (ZDGs) coupled with event-shape engineering, we identify the selection 
purity of body-body and tip-tip events in the full-overlap U-l-U collisions. With additional constraints on 
the asymmetry of the ZDG signals one can further increases the probability of selecting tip-tip events in 
U-kU collisions. 


1 Introduction 

High energy collisions between heavy ions are used to probe emergent phenomena in Quantum Chromo¬ 
dynamics (QCD), the theory of strong interaction. One feature of QCD is the transition from hadronic 
matter to a color-deconfined quark-gluon plasma (QGP) 0i as the temperature is increased. This transi¬ 
tion can occur in heavy-ion collisions of sufficient energy for the system to melt into a hot dense fireball of 
asymptotically free quarks and gluons. 

Relativistic hydrodynamic models have been successful in describing the dynamical evolution of QGP j^. 
Motivated as a testing ground for these models, a U+U collisions program was recommended in order to 


1 


study the unique collision geometry resulting from the prolate deformation of the uranium nucleus |5}|lT]. 
Such a program was carried out in 2012 at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
Lab (^. 

To understand the attraction of uranium, consider that the initial temperature distribution of each QGP 
droplet is controlled by two main factors; deterministic collision geometry (i.e. the shape of the overlap 
region between two nuclei), and quantum mechanical fluctuations in the nucleon positions. For spherical 
nuclei, the collision geometry is entirely a function of the impact parameter. However, in prolate deformed 
uranium, the geometry of the initial temperature distribution also depends on the relative spatial orientation 
of the two nuclei which can be described by the Euler angles between their long major axis. 

We focus in this paper on two limiting cases for fully overlapping uranium collisions. In one extreme 
we have “tip-tip” events, defined when the major axes of both nuclei are parallel to the beam direction. 
The opposite limit are “body-body” events, where the major axes of both nuclei are perpendicular to the 
beam direction and parallel to each other. We are interested in answering the question how, and with what 
precision, we can distinguish experimentally between these configurations. Their conceptual importance is 
explained in 0 . 


2 The model 


To model the initial energy density distribution of U+U collisions we employ the two-component (wounded 
nucleon/binary collision) Monte-Carlo Glauber model. We use the deformed Woods-Saxon distribution 


pir,0,(p) = 


Po 


1 + 


( 1 ) 


to sample the positions of nucleons inside a uranium nucleus. In Eq. 0 '. the surface diffusiveness parameter 
is c? = 0.44 fm and the saturation density parameter is po = 0.1660 fm“^ jw 15 . The spatial configuration 
of a uranium nucleus is deformed; we model its radius as fiS 
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where the average radius tq = 6.86 fm is adjusted in such a way that, after folding Eq. 0 with the 
finite charge radius of an individual nucleon, the resulting nuclear charge density distribution agrees with 
experimental constraints [^. We assume the uranium nucleus is azimuthally symmetric and choose the 
non-vanishing deformation parameters f32o = 0.28 and f34o = 0.093 for the quadrupole and hexadecupole 
deformations along its main axis, respectively. The choices of these parameters agree well with a recent 
reanalysis in [^, except for / 32 o for which Ref. [T^ gives the value 0.265. 

We use the Woods-Saxon density (Q to Monte-Carlo sample the nucleon centers and represent each 
nucleon in the transverse plane by a gaussian areal density distribution about its center: 


Pn{r±) 


^ p-rV(2B) 
(27rB)3/2 


(3) 


The width parameter B = cr)^jy(y^sj^)/14.30 depends on collision energy as described in p^ . The sum 
of these gaussian nucleon density distributions represents the nuclear density distribution for the sampled 
nucleus at the time of impact and is used to compute the initial energy density distribution generated in the 
collision. For this calculation, the two-component Monte-Carlo Glauber model weighs a relative contribution 
from binary collisions Nb and wounded nucleon participants Np j^. 

The binary collision term counts the entropy deposited by pairs of colliding nucleons and is modeled by 
a gaussian distribution with the same size as a nucleon (see Eq ■ 0 ) (^ ; the total binary collision density 
per unit transverse area is 

nBc{r±) = ( 4 ) 

where the sum is over all pairs of colliding nucleons and the normalization 7 ^ ^ is a T-distributed random 
variable with unit mean that accounts for multiplicity fluctuations in individual nucleon-nucleon collisions. 
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Each struck nucleon is said to be wounded by (or participating in) the collision and contributes a portion 
of the initial entropy density distributed symmetrically about its center; the resulting total wounded nucleon 
density per unit area is given by 

nwN(?±) = (5) 


where the sum is over all wounded nucleons in both nuclei and 7 ^ is again a fluctuating normalization factor 
with unit mean. 

We model multiplicity fluctuations in a single nucleon-nucleon collision by taking the normalizations 7 ^ j 
and 7 i to be F-distributed random variables with unit mean and variances controlled by parameters 0-qg and 
0WN. respectively. The generic F distribution with unit mean and scale parameter 0 is given by 


r(i/D)DV» - 


i/e-ip-7/e 


The multiplicity fluctuations from wounded nucleons and binary collisions are related by requiring 

1 — a, 


^pp 


-0WN = cr^BC- 


( 6 ) 

(7) 


where the parameter 0pp = 0.9175 has been fit to the measured multiplicity distributions in p-hp collisions 



The distribution in the transverse plane of the deposited entropy per unit volume is determined by mixing 
the binary collision and wounded nucleon sources using 


so(r±) 


Ks 

To 



fT-WN(f_L) + a«Bc(f_L) 


( 8 ) 


where tq is the starting time for the (hydro)dynamical evolution of the collision fireball. We choose Ks = 17.16 
and the mixing ratio a = 0.12 to reproduce the measured charged multiplicities and their dependence on 
collision centrality in Au+Au collisions at 200 AGeV. The shape of the resulting energy density distribution 
in the transverse plane is calculated from the entropy density using the equation of state (EoS) s95p-vO-PCE 
from Lattice QCD |^ . The initial energy profile is evolved using the viscous relativistic fluid dynamic code 
package iEBE-VISHNU with specific shear viscosity 77 /s = 0.08. Simulations begin at time tq = 0.6fm/c 
and decouple at a temperature Tdec = 120 MeV. The single particle momentum distribution is then computed 
using the Cooper-Eyre Formula. A full calculation of charged hadron observables that includes all hadronic 
resonance decay processes on an event-by-event basis is numerically costly; for this reason we computed 
only the directly emitted positively charged “thermal pions”, 7 r+, and take this quantity as a measure for 
total charged multiplicity. At a fixed freeze-out temperature of 120 MeV, the two quantities are related by 
a constant factor 4.6, dNch/drj ~ 4.6 (iA/r+ /dy. 

The initial energy density profiles fluctuate from event to event. Each profile can be characterized by the 
r"-weighted eccentricity coefFicients e„ and their associated "participant plane angles” <i)„: 


/ dfj_r"e“^e(r_i) 
J dr±r"e(r±) 


where (r,<p) are the standard polar coordinates in the transverse plane and e(fj^) is the initial energy density 
(^[^ . Through the hydrodynamic evolution, these spatial eccentricities {e„, 4>„} translate themselves into 
the anisotropic flow coefFicients {vn.'^n} |24] - |Z^ : 

V . .. J dippdpTe^^'^^dN/{pTdpTd(pp) 

/ d(ppdpTdN/(pTdprdpp) 


Apart from the Monte-Carlo Glauber model, there exist various other initialization models. These include 
the IP-Glasma model |^ , the MC-KLN model ||^[30] , and the TRENTO model j^. As we will see, U+U 
collisions can provide experimental measurements to distinguish between these various initializations. 
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Figure 1: Panels (a,b) show the event-averaged eccentricities £„, before hydrodynamic evolution, panels 
(c,d) the event-averaged flows Vn after hydrodynamic evolution. The left panels (a,c) represent 35,000 
minimum bias events that include multiplicity fluctuations whereas the right panels (b,d) were obtained 
from a different set of 35,000 multiplicity-selected events covering the 0-5% centrality range without 
multiplicity fluctuations. 

3 Constraining collision geometry with multiplicity, flow, and 
ZDC cuts 

3.1 Eccentricity and flow coefficients as a function of multiplicity 

In Fig. [T] we present the centrality dependence of the initial eccentricities and the final anisotropic flow 
coefficients of thermal pions for harmonic order n = 2 — 5 in U-fU collisions at 193 AGeV. In Figs. [T^,c 
minimum bias results are shown as functions of the thermal pion yields, dN^+/dy. We notice that the 
variance of £ 2,4 and V 2,4 in “most central” (i.e. highest multiplicity) collisions are larger than in the rest of 
the centrality range. This is because there the two uranium nuclei are colliding almost centrally (i.e. with 
impact parameter 6 « 0) but, as a result of the large spatial deformation, not always with full overlap. A 
mixture of tip-tip and body-body collisions in these high-multiplicity events increases the variance of the 
initial £24 which then drives a larger variance in 112 , 4 . 

In Figs.[^ and[^, we increase the statistics and focus on the 0-5% most central U+U collisions. We 
find a “knee” structure in the high multiplicity regime (< 0.5% centrality) for both £2 and V 2 - This can be 
understood as follows: First, while the ellipticity in the transverse plane for a tip-tip collision is small (as 
the overlap area is approximately circular), body-body collisions produce an ellipsoidally deformed overlap 
region with larger ellipticity £ 2 . Second, although fully overlapping tip-tip and body-body collisions share the 
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Figure 2: The ellipticity as a function of dS/dy from the MC-Glauber model, for collisions roughly in 
the 0-5% centrality range, with (blue dashed line) and without (black solid line) including multiplicity 
fluctuations from single p-l-p collisions. 


same number of participants, more binary collisions between nucleons can happen in the optically thicker 
tip-tip event, implying (in our two-component Glauber model) a larger initial dS/dy deposited in the tip-tip 
configuration. In the presence of fluctuations which lead to a range of £2 values for a given dS/dy and 
vice-versa, the larger average multiplicity in tip-tip collisions implies an increasing bias toward small £2 when 
selecting events with larger and larger values of dS/dy. This preferential selection of tip-tip orientations at 
high multiplicities accounts for the appearance of a knee structure in the initial ellipticity (Fig. Q. We 
see in Fig. [J that the knee is preserved after an event-by-event hydrodynamic simulation when plotting the 
elliptic flow of the final particle distribution as a function of multiplicity. 

We emphasize that experimental results from STAR do not show this knee structure (I^ . Considering 
the preservation of the structure after hydrodynamic evolution as seen in Fig.j^ we conclude that, in contrast 
to Au-hAu collisions where it has been extensively tested, the two-component MC-Glauber model fails to 
correctly identify entropy production in ultra central U-FU collisions where the knee is predicted by the model 
but not found experimentally. Hence the non-linear dependence of multiplicity on the number of wounded 
nucleons observed in spherical Au-fAu and Pb+Pb collisions as a function of collision centrality cannot be 
attributed to a binary collision component as implemented in the two-component MC-Glauber model. 

Some corrections to the entropy production in these ultra central events arise from the inclusion of p+p 
multiplicity fluctuations. We see in Fig. that adding said multiplicity fluctuations weakens but does not 
erase the knee structure in £2 vs. dS/dy. Hence, this effect alone does not appear sufficient to reach 
agreement of the MC-Glauber model with data for ultra central U+U collisions. We acknowledge that 
more drastic fluctuation models have been suggested in order to more successfully adjust the theoretical 
predictions of MC-Glauber to experimental results. 

Also of note is the success of the gluon saturation physics as implemented in the IP-Glasma model. 
Interestingly, this model is able to simultaneously accomodate a strong nonlinearity of dN/dy as a function 
of A^part in Au-FAu and Pb+Pb and a weak dependence of dN/dy on collision orientation in central U+U at 
fixed number of participants, while the MC-Glauber model cannot j^. Alternatively, it has been suggested 
in |3^ that a model that produces entropy according to the number of wounded valence gluons (rather 
than wounded nucleons) can also reproduce the observed nonlinearity of dN/dy as a function of particpant 
nucleons in Au+Au and Pb+Pb at RHIC and LHC, without a binary collision component. It would be 
interesting to study the production of such a model for central U+U collisions of varying orientations. 
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Figure 3: Probability distribution for dS/dy (left) and S2 (right) for different event classes within a 
sample of 1% ZDC events. The top panels show the distributions of tip-tip and body-body collisions 
scaled according to their contribution to the total population within the 1% ZDC sample. The bottom 
panels show relative probabilities for tip-tip and body-body events among all events of a given dS/dy 
(left) or £2 (right). 


3.2 Selecting high overlap events with combined ZDC and multiplicity cuts 

In an experimental analysis of relativistic heavy ion collisions, the charged hadron multiplicity, dN^h/dy and 
its elliptic flow coefficient V 2 can be used to classify events. Hydrodynamic studies have shown that the 
initial £2 maps linearly to the V 2 of hadrons and the the initial dS/dy is monotonically related the final 
total particle multiplicity, dN/dy [^. We can therefore use dS/dy and £2 from the initial conditions as a 
satisfactory proxy for charged hadron dN^h/dy and V 2 to test whether we can select the fully overlapping 
tip-tip and body-body U+U collisions. 

For our analysis, we make theoretical approximations for the use of experimental forward and backward 
zero degree calorimeters (ZDCs). Placed at zero degrees far from the colliding pair, ZDCs catch information 
about the spectator neutrons that pass through a collision without participating. We classify our collisions 
by using the number of spectators Ns = 476 — A^part to mimic the experimental ZDC signal j^. For our 
study we look at 65,000 events in the 1% most participating ZDC range {Ng < 19). Selecting on the most 
participating ZDC collisions allows for a restriction of the set of collisions to more fully overlapping events. 
In such a regime, any initial geometric effects should come more exclusively from the deformed shape of the 
uranium nucleus. 

We define the tip-tip and body-body event classes using the pair of angles (di, 2 , (/i’l.z) from the two 
incoming nuclei, where 9 denote the polar angle between the long major axis of the uranium nucleus and 
the beam direction and (f) is the azimuthal angle in the transverse plane. An event is defined as tip-tip if 
>/cos^ 01 -F cos^ 6*2 > 0.86. We classify an event as a body-body event if both Vcos^ 6*1 -F cos^ 6*2 < .31 and 
\4'i ~ '/' 2 I < tt/IO. The polar angle constraints imply that for equal 61 = 02 , this common angle 0 is less 
than tt/IO for tip-tip and greater than 47 r /10 for body-body. For body-body events, the additional azimuthal 
constraint forces alignment of the long major axes. 

In Fig. I^we plot the probability distributions for £2 and dS/dy. Using our collision definitions we can 
directly read off from the figure the likelihood of selecting a certain orientation based on a given eccentricity 
or multiplicity cut. We see in the left bottom panel that by cutting (within our 1% ZDC sample) on events 
with large dS/dy we can enrich the fraction of tip-tip events to about 50%, whereas cutting on low dS/dy 
enriches the fraction of body-body events, but never to more than about 20%. The 20% limit arises from 
admixtures from imperfectly aligned collisions that are not really ’’full overlap”. The enrichment of tip-tip or 
body-body by varying dS/dy relies on the assumed two-component nature of entropy production which also 
produced the knee structure discussed before. Indeed, selection efficiency of specific collision geometries by 
cutting on dS/dy is model dependent. 

We therefore consider “event engineering”, i.e. selecting events by the magnitude of their V 2 flow vectors 
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Figure 4: The black curve is the distribution of tip-tip events scaled according to their contribution to 
the total probability distribution for £2 as seen in the bottom right pannel in Fig. The blue dashed 
curve shows the increased contribution of tip-tip collisions within the 10% of events having the smallest 
difference in participants ATVpart (a proxy for ZDC correlation). 


(for us, of the linearly related £ 2 ), shown in the right panels of Fig. Since tip-tip events have on average 
smaller ellipticities (see upper right panel), selecting events with small ellipticity (or, in experiment, small 
V 2 ) enriches the tip-tip fraction. However, in this way we will never reach more than about 25% purity of 
the tip-tip sample. On the other hand, cutting the 1% ZDC events on large £2 (or V 2 ) will enrich the sample 
in body-body events, with a purity that can reach about 40% for the largest £2 values. While we have not 
yet been able to verify this with an actual cut on V 2 , we expect this feature to survive the hydrodynamic 
evolution due to the almost perfect linearity between £2 and V 2 - 

The current ZDC cut strategy can be refined further to increase the probability of selecting tip-tip events. 
Rather than looking at the ZDC signal in one of the two ZDC detectors or the sum of the ZDC signals in 
both detectors, we can look at the correlation of these two signals. Events with equal forward and backward 
ZDC signals (i.e. equal numbers of spectators from both nuclei) provide a better definition of the categories 
full overlap, tip-tip, and body-body than events with asymmetric ZDC signals where all spectators come 
from only one of the colliding nuclei. The difference in participants AApart = |.^part,i — .5i’part.2| quantifies 
the ZDC correlation in our model. Low values of ANp correspond to the most correlated forward and 
backward ZDC signals. To demonstrate one application, we reconsider the events from the bottom right 
panel in Fig. |^and now select from the sample only events in the lowest 10% of AApart- The selection 
on small values of AlVpart eliminates from the sample asymmetric configurations that we loosely describe 
as “tip-body”. Collisions of this type produce low values of £2 without the angular criteria necessary to 
be considered tip-tip and therefore dilute the contribution of the true tip-tip configurations at the lower 
range of £ 2 . We show in Fig. |^that selecting the lowest 10% of AApart increases the selection efficiency 
of £2 for tip-tip configurations by a factor of 1.4. As a final comment, we point out that it might also be 
interesting to use ZDC correlations in the opposite way and to select and study events with asymmetric 
tip-body configurations. 
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4 Conclusion 


Within the two-component MC-Glauber model for initial energy production, the prolate deformation of the 
uranium nucleus was shown to generate a knee in the centrality dependence of the ellipticity of the initial 
temperature distribution. The knee was seen to be preserved by hydrodynamic evolution, after which it 
manifests itself in the centrality dependence of V 2 - Such a knee structure is not seen in the STAR data. This 
rules out the two-component MC-Glauber model for initial energy and entropy production. An enrichment 
of tip-tip configurations by triggering only on high-multiplicity in the U+U collisions thus does not work. 

To increase the selection capability between different collision geometries, we impose combined cuts 
on initial conditions using the spectators (ZDC), dS/dy, and 62 - For 1% ZDC events, we found that we 
could enrich tip-tip collision geometries to about 50% by cutting on high multiplicity within that sample, 
and body-body configurations to about 20% purity by selecting low-multiplicity events. These numbers rely 
on the binary collision admixture in the two-component MC-Glauber model and are thus model-dependent. 
They do include effects from multiplicity fluctuations. 

We also studied the efFiciency of selecting different collision geometries by “event engineering”, i.e. by 
cutting on 82 (by cutting on V 2 in the experiment). In this case events selected for high 62 can enrich 
body-body collisions to about 40% purity while cutting on low £2 gives a tip-tip sample with about 25% 
purity. The latter can be boosted to about 35% purity by eliminating events with asymmetric ZDC signals. 
These results should not be sensitive to the binary collision admixture in the two-component MC-Glauber 
model and thus should be less model dependent. 
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